Eigenwavelets of the Wave Equation 



Gerald Kaiser 

Signals & Waves, Austin, TX • www.wavelets.com • kaiser(3wavelets . com 



Summary. We study a class of localized solutions of the wave equation, called 
eigenwaveJets, obtained by extending its fundamental solutions to complex space- 
time in the sense of hyp erfunct ions. The imaginary spacetime variables y, which form 
a timelike vector, act as scaJe parameters generalizing the scale variable of wavelets 
in one dimension. They determine the shape of the wavelets in spacetime, making 
them pulsed beams that can be focused as tightly as desired around a single ray by 
letting y approach the light cone. Furthermore, the absence of any sidelobes makes 
them especially attractive for communications, remote sensing and other applica- 
tions using acoustic waves. (A similar set of 'electromagnetic eigenwavelets' exists 
for Ma:xweirs equations.) I review the basic ideas in Minkowski space R^'^, then 
compute sources whose realization should make it possible to radiate and absorb 
such wavelets. This motivates an extension of Ifuygens' principle allowing equiva- 
lent sources to be represented on shells instead of surfaces surrounding a bounded 
source. 



1 Extension of wave functions to complex spacetime 

The ideas to be presented here affirm tfiat complex analysis resonates deeply 
in "real" physical and geometric settings, and so they are close in spirit to 
the work of Carlos Berenstein (see BG91 , BG95, B98 for example), to whom 
this volume is dedicated. 

Acoustic and electromagnetic wavelets were first constructed in |K _94| . It 
was shown that solutions of homogeneous {i.e., sourceless) scalar and vector 
wave equations in Minkowski space R^'^ extend naturally to complex space- 
time, and the wavelets were defined as the Riesz duals of evaluation maps 
acting on spaces of such holomorphic solutions. The sourceless wavelets then 
split naturally into retarded and advanced parts emitted and absorbed, re- 
spectively, by sources located on branch cuts needed to make these parts 
single- valued. Later work |K31 IK4| was aimed at the construction of realizable 
source distributions which, when synthesized, would act as antennas radiating 
and receiving the wavelets. Two difficulties with this approach have been (a) 
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that the computed sources are quite singular, consisting of multiple surface 
layers that may be difficult to realize in practice, and (b) in the electromag- 
netic case the sources appeared to require a nonvanishing magnetic charge 
distribution, which cannot be realized as no magnetic monopoles have been 
observed in Nature. In this paper we resolve the first difficulty by replacing the 
spheroidal surface supporting the sources in (KlEll by a spheroidal shell. It 
is shown in K4a that the second difhculty can be overcome using Hertz po- 
tentials, which give a charge-current distribution due solely to bound electric 
charges confined to the shell. 

Although our constructions generalize to other dimensions, we shall con- 
centrate here on the physical case of the Minkowski space W^'^. Let 

X = ir,t), y = {a,b) e R""-^ (1) 
be real spacetimc vectors and define the complex causal tube 

T = {x- iy eC^ : y is timelike, i.e., \b\ > \a\}. (2) 
It was shown in |K94[ IK3j that solutions of the homogeneous wave equation 

nfo{x) = {d?-A)Mr,t) = (3) 
extend naturally to analytic functions fo{x — iy) in T in the sense that 

lim \foix- iy) - foix + iy) \ = foix), (4) 

where y —^ +0 means that y approaches the origin within the future cone, 
i.e., with b > \a\. This kind of extension to complex domains is familiar in 
hyperfunction theory; see jK88[ IKS99j for example. We now show that even 
when the wave function has a source, i.e., 

Dfix)=A7rg{x), (5) 

it extends analytically to T outside a spacetime region determined by the 
source. It will suffice to do this for the retarded propagator 

Sjt - r) 

G{x) = , (6) 

r 

which is the unique causal fundamental solution: 

□G(.x) = A7rS{t)6{r) = AnS{x), G{r, t) = Vt < 0. (7) 

If the source g is supported in a compact spacetime region VF, the unique 
causal solution of Q is given by 



fix) = / dx' G{x~ x')g{x'). 
Jw 



(8) 
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Assume for the moment that G{x) has been extended to G{x — iy). Then we 
define the source of G as the distribution S in real spacetime given by 

4TrS{x - iy) = U.^G(x - iy), (9) 

where D^; means that the wave operator acts only on a;, in a distributional 
sense, so that the imaginary spacetime vector y is regarded as an auxiliary 
parameter. The extended solution is now defined as 

j{x - ^y) = / dx' G[x - x' - iy)g{x') (10) 
Jw 

and it satisfies the wave equation 

^xRx - iy) = iTTg{x - iy) 

with the extended source 

g{x-iy)^ dx' 6{x ~ x' - iy)g{x'). (11) 
Jw 

Formally, the extended delta function 5{x — iy) is a 'point source' at the 
imaginary spacetime point iy as seen by a real observer at x. Actually, it will be 
seen to be a distribution in x with compact spatial (but not temporal) support 
localized around the spatial origin (r — 0) and depending on the choice of a 
branch cut needed to make G single-valued. This branch cut is precisely the 
region where G fails to be analytic, and the integral H10|) determines a region 
W containing W where / fails to be analytic. 

A general solution fi{x) of © is obtained by adding a sourceless wave 
/o(x) to (O. Since /o is analytic in T, /i is analytic in T outside of W. It 
therefore suffices to concentrate on the propagators as claimed. In the rest 
of the paper we construct extended propagators, study their properties, and 
compute their sources. 

2 Extended propagators 

In accordance with (QJ, we use the following notation for complex space and 
time variables: 

r = r-meC^, i = t - ib e C 
X ^ X ~ iy ^ {f,i) e T 4^ |6| > |a|. 

As above, we interpret ia formally as an imaginary spatial source point, so 
that f is the vector from the imaginary source point ia to a real observer at 
r. To extend the propagator ((HJ, begin by replacing the one-dimensional delta 
function with the Cauchy kernel, 
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6{t) ^ ~S{i) ^ t^t-ib, (12) 
2mt 



which indeed satisfies a condition of type Q : 

^hni^ ^5{t - ib) - ~5{t + = 5{t). (13) 



To complete the extension of G{r,t), we must also extend the Euchdean dis- 
tance r(r) — \r\. Define the complex distance from the source to the observer 
as 



f(r) = Vf~f = \/r'^ — o? — 2ir ■ a, where r = |r|, a = |a|. (14) 

f(r) is an analytic continuation to C'^ of r(r). Being a complex square root, 
it has branch points wherever f • f = 0. For fixed a 7^ 0, these form a circle 
of radius a in the plane orthogonal to a,^ 

C = {r eM.^ : r = 0} = {r : r = a, r ■ a^Q}. (15) 

To be consistent with the notation f = r ia, we write 

r = p — iq. (16) 

Comparison with (|14|l gives the following relations between {p, q) and the 
spherical and cylindrical coordinates with axis along a: 

p^ — — — a^, pq — a ■ r = ar cos 6 = az (17) 

and 

a^p^ = a\r^ - z') - a\a' + / - q') - p\' 

^{a'+p')(a'-q'). (18) 

It follows that the real and imaginary parts of r are bounded by r and a, 
respectively: 

< r^, i.e., |Ref|<|Ref| 
q'^<a^, i.e., |Imf|<|Imf|, (19) 

with equalities attained only when r is parallel or antiparallel to a. 

Since a will be a fixed nonzero vector throughout, we will usually regard 
f,p,q as functions of r only, suppressing the dependence on a. Note that 
R'^ — C is multiply connected since a closed loop that threads C cannot be 
shrunk continuously to a point without intersecting C. In particular, if we 
continue f analytically around a simple closed loop, we obtain the value — f 
instead of f upon returning to the starting point. Thus f is a double- valued 

^ In R", C would be a sphere of codimension 2 orthogonal to a. 
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function on M^. To make it single- valued, we choose a branch cut that must 
be crossed to close the loop. Instead of returning to the starting point as — f, 
the sign reversal now takes place upon crossing the cut. To give an extension 
of the positive distance, the branch must be chosen so that 

a^O f ^ +r, (20) 

and the simplest such choice is obtained by requiring 

Re f = p > 0. (21) 

The resulting branch cut consists of the disk spanning the circle C, 

V = {r£R^:p = Q} = {r: r <a, r ■ a^O}, dV^C. (22) 

T) will be called the standard branch cut and f the standard complex distance. 
General branch cuts, obtained by deforming V while leaving its boundary 
intact, will be considered in the next section. 

If the observer is far from C, it follows from (|14|l and H2U|I that 

r^a p Ki r and q^acosO, (23) 

Thus, {p, q/a) are deformations of the spherical coordinates (r, cos 9) near the 
source. From (|17|) and (|18|l it follows that level surfaces of p^ (as a function 
of r, keeping a fixed) are spheroids Sp and those of arc the orthogonal 
hyperboloids Hq, given by 

2 2 

^P- ^^ + ^ = 1' (24) 
+ p"^ 

Hq : ' , --7 = 1, 0<q^ <a^. (25) 
q^ — q'^ 

All these quadrics are confocal with C as the common focal set. As p — > 0, 
Sp collapses to a double cover of the disk V. The variables {p,q), together 
with the azimuthal angle 4> about the a-axis, determine an oblate spheroidal 
coordinate system, as depicted in Figure 

We now define the extended propagator as 

G(f,0 = M^ = --^. (26) 
r 2iTir{t — r) 

This is our basic wavelet,'^ from which the entire wavelet family is obtained 
by spacetime translations: 

G^{x)^G{x- z), z = x' + iy={r' + ia,t' + ib) eT (27) 



^ In applications, it is better to use time derivatives of G, which have vanishing 
moments and better temporal decay and propagation properties |K4| . 
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The family Gz may be called eigenwavelets of the wave equation in the 
sense that they are proper to that equation, though of course they are not 
eigenfunctions. In fact, Gz{x) is seen |K3| to be a pulsed beam originating 
from r — r' at t — t' and propagating along the direction of a/b, i.e., along 
a if y is in the future cone and along —a if y is in the past cone. The pulse 
has a duration |6| — a along the beam axis. By letting y approach the light 
cone {a —^ the beam can be focused as tightly as desired around its axis, 
approximating a single ray along y. Equation (|10|l states that the extended 
causal solution f{x — iy) is a superposition of eigenwavelets, all with the same 
y. This gives a directional scale analysis of the original solution /(x) which 
may be called its eigenwavelet transform. 

The eigenwavelets have the spheroids Sp as wave fronts and propagate out 
along the orthogonal hyperboloids Ti.q with strength decaying monotonically 
away from the front beam axis. Hence they have no sidelobes, which makes 
them potentially useful for applications to communication, radar and related 
areas. These properties are illustrated in Figures [21 and 13 

We may visualize the effects of the extension G(r, t) — > G{f, t) as follows. 
The extension t ^ t replaces the spherical impulse S{t — r) in by a spherical 
pulse S{t^r) of duration \b\. The extension r ^ f then deforms this spherical 
pulse to a pulsed beam in the direction of a/b. By H23|l. 

r ^ a ^ f Ki r — ia cos 9, (28) 

hence the larger we choose a, the stronger the dependence of f on cos 9 in the 
far zone and the more focused the beam. 

Let us emphasize that G depends on the complex spatial vector r e 
only through the complex distance f by writing 



W{f,t) = G{f,t)^——^—. 

ZTnr[t — r) 



(29) 
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Fig. 2. Time-lapse plots of \G{x — iy)\ in the far zone, showing the evolution of 
a single pulse with propagation vector y = (0,0, 0,6). Clockwise from upper left: 
b/a — 1.5, 1.1, 1.01, 1.0001. As b/a 1, y approaches the light cone and the 
pulsed beam becomes more and more focused around the ray y. We have taken the 
slice X2 ~ 0, so that the disk "D becomes the interval [—a, a] on the Xi-axis and the 
pulse propgates in the direction of the xi-x^ plane. 

Due to the factor f in the denominator, if" is discontinuous across V and 
singular on C. T) generalizes the point singularity of G at r = and will be 
the spatial support of the source l^. To avoid further singularities, the factor 

t~f= [t — p) — i{h — q) 

must not vanish for any r. By (|19|l . 

b-q^Q^r ^ a< (30) 

so a necessary and sufficient condition for !f'(r, 1!;) to be analytic whenever 
r ^ V is that (f , t) G T. Recalling that the tightness of the beam is controlled 
by the size of a, ()30|l means that the beam cannot become tighter than a 
single ray and, in fact, fails to be analytic along the ray in the limit a = \b\. 
The volume element in R'^ in oblate spheroidal coordinates is 

dV = -{p^ + q^)dpdqd(p = - \f\^dpdqd(P, (31) 
a a 

hence 'F is locally integrable and square integrable. A differentiation gives 

AttS{x - iy) EE a^G{x -iy) = (x -iyeT, r ^ V). (32) 

Therefore Slx — iy), with y a fixed timelike vector, is a distribution in x — (r, t) 
with spatial support in V. (The temporal support is noncompact; in fact, 
6{x — iy) decays as 1/i due to the Cauchy kernel.) 
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Fig. 3. Near-zone grapiis of \G{x — witli y = (0,0, 1, 1.01) immediately after 
launch, evolving in the xi-xs plane with 0:2 = as in Fig. 1. Clockwise from upper 
left: t = 0.1, 1, 2, 3. The ellipsoidal wave fronts and hyperbolic flow lines are visible. 
The top of the peak is cut off to show the behavior near the base. The two spikes 
represent the branch circle, whose slice with X2 = consists of the points (±1, 0, 0). 

The source 5{x — iy) was computed explicitly in |K3| and turns out to 
be quite singular. It consists of a single layer and a double layer on 2?, both 
of which diverge on the boundary C where If^ is singular. We will compute 
regularized versions of W and 5 by using the freedom to deform the branch 
cut to eliminate the singularity on C. 

3 Regularization by branch cut deformation 

A general branch cut ;S is a membrane obtained by a continuous deformation 
of the disk V leaving its boundary intact, 

dB^C. (33) 

B inherits an orientation from V, which in turn is oriented by a. Let Vg be 
the compact volume swept out in the deformation from V to B. Let us define 
the complex distance with branch cut B in terms of f = r^, by 
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r. = r (34) 
I — f if r e Ve • 

I claim that fg is continuous except for a sign reversal across B, generalizing 
the sign reversal of f across V. This can be seen most simply if B does not 
intersect the interior of T), so that they have only the boundary in common. 
Then Vs is either all on the 'positive' or all on the 'negative' side of V. If Vg 
is 'positive,' then its boundary is 

dVB=B- V, (35) 

meaning that the orientation (outward normal) of the boundary is positive 
on B and negative on T). Since f changes sign upon crossing V 'upward' 
into Vb, I|34|I shows that fg is continuous across "D. This proves that its only 
discontinuity is the sign reversal in crossing B, as claimed. Similarly, if Vb is 
'negative,' then its boundary is 

dVB=V-B (36) 

and the above argument remains valid. To handle branch cuts that intersect 
the interior of I?, we restate the 'negative' case (j3fi|l by declaring Vb negatively- 
oriented, so that its boundary is oriented by the inward normal. Denoting the 
negatively oriented volume by —Vb, H36|) can be restated as 

d{-VB)=B-V. (37) 

Hence the rule H35|l applies to every branch cut B obtained by a continuous 
deformation of V, whether or not it intersects the interior of V, provided the 
orientation of the swept-out volume is taken into account. If B intersects the 
interior of V, then Vg has both positively and negatively oriented components. 
The definition (|34() of the branch fg remains valid whether Vb (or any of its 
components) is positively or negatively oriented. 

Of special interest will be the upper and lower spheroidal branch cuts 

S ± = U Aa, (38) 

where denote the upper and lower hemispheroids 

5± = {r e 5„ : ±z > 0} 

and 

Aa= {r : r ■ a = 0, < r'^ < a'^ + a^} 

is the apron connecting them to C, which must be included so that dB^ = C 
as required. The cut S+ is depicted in Figure^ 

We can now construct a regularized version of the extended propagator 
by taking the average of the propagators with cuts B^ and B^ . Denote the 
complex distances with cuts B^ by r± instead of fg ± , and let 
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Fig. 4. Tlie upper liemispiieroidal brancli cut witii its apron. 



tf'^(r~,i) = i =GAix~ ly). (39) 

Let be the interiors of the upper and lower hemispheroids. By (|34ll . 

r ^ f + = -f, f_=f (40) 

T" € =4> f+ = f, f_ = — f. (41) 

Hence, in both Vj" we have 

, f) = ^ I = J- , (42) 

All if it — f) 47rif(i + f) 27ri(t^— f^) 

which is independent of the choice of branch. This shows that the disconti- 
nuities across the aprons cancel in the average 'I'a- Furthermore, by (|30|l we 
have 

|6| > a P - = (t- f)(t + f) 7^ 0, 

showing that the singularities on C cancel as well. That is, is analytic at 
all interior points of the spheroid Sa ■ 

In the exterior of Sa we have r± — f and hence <Fa = Since T) is 
contained in Sa and 'F is analytic outside of T), we conclude that 'FA{r, i) fails 
to be analytic only when r E Sa- Denoting the interior field by 'Fi and the 
exterior field by we have 

tf'i (f , i) = i i) + tf'(-r~, i) } (43) 

^2{f,i) = >f(f,t). 

Thus ^A is analytic except for a bounded jump discontinuity 3,cross Set §iv6ii 

by 

tf'j(f, i) = i) - <fi(f , i) = i {<f'(7^, i) - ^{-f, i)} = Gj(x - ly). (44) 
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It follows by the same argument as in (|32|) that the source distribution 

4TrSAix-iy) = n^^A{f,i) (45) 

is supported spatially on Sa ■ Because Sa is obtained by twice differentiating a 
discontinuous function, it consists of a combination of single and double layers 
on Sa- But the jump discontinuity in 'Fa is bounded (unlike that in tZ^, which 
diverges on C), and so are these layers; see |K4| . 

The above arguments remain valid if instead of we use any two branch 
cuts whose common interior V contains the branch circle C. In that case, the 
averaged propagator is analytic in T except for a finite discontinuity when r 
crosses the boundary dV, and its source distribution is supported spatially 
on dV. However, the above choice has the advantage that dV = Sa are wave 
fronts, hence all parts of the surface radiate simultaneously and coherently. 

4 Extended Huygens sources 

Let H be the Heaviside step function. Since < p < a in the interior of Sa 
and p > a in the exterior, we have 

i'Air, i) = H{a - p)1'i{f, i) + H{p - a)W2{f, i) (46) 

where the interior and exterior fields are given by 143(1 . This can be used to 
compute the source distribution 5a defined in 145|l . and the result is sum of 
terms with factors 5{j) — a) and 6'{p — a). The former are interpreted as single 
layers on Sa, and the latter as double layers. 

An interesting practical question is whether the wavelets Wa, interpreted 
as acoustic pulsed beams, can be realized by manufacturing their sources. A 
similar question can be posed for their electromagnetic counterparts, which 
solve Maxwell's equations; see K4 . It is doubtful whether an acoustic source 
can be produced including double layers, and the problem becomes even more 
difficult in the electromagnetic case because the current density involves yet 
another derivative, hence a still higher layer |K4a| . The multilayered structure 
is unavoidable as long as we insist on surface sources. We now propose a 
method for constructing solutions of the wave equation where the transition 
occurs in a shell instead of a surface. It will be simpler to present this method 
initially in a somewhat more general context. 

Given a function p(r, t) on R"'"'^ and two regular values pi < p2 in its 
range, define two time-dependent surfaces and volumes in M" by 

Si{t) = {r : p{r, t) = pi}, S2{t) - {r : p(r, t) = p^} 
Vi{t) - {r :p(r,t) <pi}, V2{t) - {r : p{r,t) > P2}. 

Let /i, /2 be solutions of the wave equation in R"'^ with sources 51,32: 
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□/fc(r,t)=gfe, k^l,2. (47) 
We want to construct an interpolated solution ,f(r,t) such that 

f{r,t) = fkir,t)yreVk{t) (48) 
and compute its source. This can be done by choosing functions hk{r, t) with 

h,{r,t) = \]; ''l^'^-^^, h,{r,t)^l^h,{r,t) (49) 
[0, reV2{t) 

and letting 

/ - hih + /i2/2 = hkfk (50) 

where the (Einstein) summation convention is used. The source of / is found 
to consists of two parts, 

g = Df^g,+g^, (51) 

where 

gi = hkgk (52) 

is an interpolated source and 

g^ = 2hkh - 2Vhk ■ Vfk + iDhk)fk if = dtf) (53) 

is a transitional source which, by H49|l . is supported on the transition shell 

VT{t) = {r:pi<p{r,t)<p2} (54) 

and depends only on the jump Geld f , = /2 — ,fi ■ 

gr = 2h2fj - 2V/i2 • V/,, + {nh2)fj. (55) 

Now suppose that Vi(t) and Vrit) are compact and we are given only one 
source g2, supported in Vi(t). Letting /2 be its causal field, our objective 
is to find an equivalent source g supported in Vrit) whose causal field / is 
identical with /2 in V2 (t) . It suffices to choose any solution /i whose source 
gi is supported in V2{t), since the interpolated source (|52(l then vanishes and 
hence g ^ gr- fi is a sourceless internal Geld in Vi{t), and the source so 
constructed on Vrit) generalizes the idea of a Huygens source on a surface 
surrounding the support of 52 ■ We may recover the latter by assuming that p 
is time- independent (hence so are Sk and Vk) and choosing hk{r) so that 

lim V/i2('") — '5(p(r) ~ p2)n{r) 
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where n{r) is a field of orthogonal vectors on 5*2 pointing into V2. The corre- 
sponding scheme in the electromagnetic case reduces to the usual boundary 
conditions on an interface between two media |K4a] . 

Returning to n = 3 with p — Re f, let fk = as in (|43|l and hk be 
time-independent {e.g., functions of p only). A smoothed version of 'Fa H39() 
is 

tZ/™ = hi<Pi + h2l'2- (56) 
Since 'Fk are sourceless in Vr, (|51|l gives the smoothed version of H45|l as 

47r5r = O^'^T = -2V/i2 • V^'j - {Ah2)'Pj (57) 

where 

is the jump field from Vi to V2 as in H44() . but no longer restricted to a single 
spheroid Sa ■ If we now let pi P2 — ct and 

hi^H{a~p), h2 = H{p-a), 

then the transition becomes abrupt on Sa and If^J™ becomes 'Pa I|39|I . Since 

V/i2 (5(p - a)Vp 

Ah2 = 6'{p - a)|Vpp + S{p - a)Ap, 

equation (|S7|l becomes 

AttSa = -2(5(p - a)Vp • V'f'j - S{p - a)Ap'Pj - S'{p - a)\\'p\'^'I'j 

displaying the aforementioned single and double layer structure on Sa • To get 
an explicit expression, use |K4I Appendix] 

pr + qa 2p 



= 2 I 2 ' "^P 2 , 2 
and 

Vp • V!f J = 'PjVp ■ Vf = iP'jiVpp 

where 'Pj is the complex derivative of ']/{f,t) with respect to f (keeping in 
mind that !P'(±f, t) are analytic in f for p > 0), 
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5 Conclusions 

Although I have concentrated on the wave equation in four-dimensional 
Minkowski space M^'^, similar considerations apply in M"'^. In fact, the awk- 
ward extension of the propagator, using the Cauchy kernel in time but the 
complex distance in space, becomes much more natural when G(r, t) is viewed 
as the retarded part of the analytic continuation of the fundamental solution 
Ge{R) of Laplace's equation in Euclidean M"^^ |KOI IK3| . based on the com- 
plex distance 

R= Vr R, Re C"+\ 

whose branch points form a sphere S"""^ in R"^^ of codimension 2 and ra- 
dius I Im ^|. The extended delta function ^^(H),'^ defined by applying the 
Laplacian in R to the extension Ge{R), is supported on 5*""^ for odd n > 3, 
but a branch cut, consisting of a 'membrane' bounded by S"~^, is needed in 
all other cases.'' Given any test function / in R""*"^, the convolution 

f{R)^ f L{R-R')f{R')dV{R') (58) 

defines an extension of / to C"^^, non-holomorphic in general, whose restric- 
tion to the Minkowski subspace R"'^, obtained by letting R = {r,it), is a 
solution of the following initial-value problem for the wave equation: 

{df-Ar)f{r,it) = (59) 
/>,0) = /(r,0) (60) 
(dt - idi)f{r, b + It) \i,^t=o= 0. (61) 

For odd n > 3, the proof of ifS^ is based on the fact that 5e is distributed 
uniformly on 5*""^ and hence / is a spherical mean of / |J55| . This relates 
the support of Se for odd n > 3 to Huygens principle. The other cases can be 
treated by applying a distributional version of Hadamard's method of descent. 

Equation (|?)T|l states that /(r, b + it) satisfies the Cauchy- Riemann equa- 
tion in its last variable; but since this holds only at one point, it does not 
imply analyticity — as it cannot since /(r, 6) need not have any analytic con- 
tinuation in b. If one exists, it is indeed given by /(r, b + it). This generalizes 
an old theorem by Paul Garabedian |(t64I pp 191-202]. 
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* This is because Ge{R) = Cn/R"~^ for n > 2 and Ge ~ ci logi? for n — 1. 
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